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The spin Hall and longitudinal conductivity of a 2D heavy-hole gas with A;-cubic Rashba and 
Dresselhaus spin-orbit interaction is studied in the ac frequency domain. Using Kubo linear-response 
theory and a recently proposed definition for the (conserved) spin current operator suitable for 
spin-3/2 holes, it is shown that the spin conductivity tensor exhibit very distinguishable features 
from those obtained with the standard definition of the spin current. This is due to a significant 
contribution of the spin-torque term arisen from the alternative definition of spin current which 
strongly affects the magnitude and the sign of the dynamic spin current. In the dc (free of disorder) 
limit, the spin Hall conductivity for only (or dominant) fc-cubic Rashba coupling is a^y (0) = — 9e/87r, 
whereas o%y{0) = — 3e/87r for only (or dominant) fc-cubic Dresselhaus coupling. Such anisotropic 
response is understood in terms of the absence of mapping the fc-cubic Rashba «-» Dresselhaus 
Hamiltonians. This asymmetry is also responsible for the non-vanishing dc spin Hall conductivity 
(crjy (0) = — 6e/87r) when the Rashba and Dresselhaus parameters have the same strength, in 
contrast with its corresponding case for electrons. These results are of relevance to validate the 
alternative definition of spin-current through measurements in the frequency domain of the spin 
accumulation and/or spin currents in 2D hole gases. 



I. INTRODUCTION 

Nowadays spintronics research is devoted to great ex- 
tent to the study of fundamental issues concerning how 
to generate, manipulate and detect spin currents in a 
controlled and efficient wayi In particular, spin currents 
can be generated in solid state systems through the spin 
Hall effect (SHE) j2j2i» a relativistic spin-orbit interaction 
(SOI) phenomenon that converts an electrical voltage 
into a pure spin current. The SHE is a result of a spin- 
dependent scattering from charged impurities through 
side jump and/or skew processes^ due to spin-orbit cou- 
pling (extrinsic SHE) 2-3 . It can be generated also in bulk 
and heterostructure due to built-in fields that modifies 
the band structure via the SOI (intrinsic SHE). 6 ' 7 

The SHE has been inferred by optical means via the 
spin Hall accumulation at the edges of a doped (n-type) 
semiconductor channel 8 ^, as well as in two-dimensional 
(p-type) hole gases 10 (2DHGs). It has been also detected 
electrically in metals throughout its reciprocal effect, the 
inverse SHE ; 11 ' 12 ! 13 - the generation of a charge current 
by a transverse spin-current in materials with SOI^ Re- 
cently evidence of a rather large spin Hall signal (2.9 mf2) 
at room temperature in FePt/Au multitcrminal devices 
exhibiting significant skew scattering processes has been 
also observed^ Most recently, Briine et alr^ have re- 
ported the first electrical measurement and manipulation 
of the intrinsic SHE using ballistic HgTe semiconductor 
nanostructures observing even larger spin Hall signals of 
the order of few fcO. 

Theoretically, the intrinsic dc spin Hall conductivity in 
2D electron systems with fc-linear Rashba and/or Dres- 
selhaus type of SOI is known to be suppressed by weak 
(non magnetic) disorder by the appropriate inclusion of 
vertex corrections. 17 However such cancelation does not 



hold in the presence of magnetic fields and/or magnetic 
impurities i 18 ' 19 Furthermore, the vertex corrections are 
in fact identically zero in 2D-hole systems with fc-cubic 
Rashba, 20 2D fc-cubic Dresselhaus electron systems^ 21 
and Luttinger Hamiltonians, 22,23 for which the spin-Hall 
conductivity is shown to be robust against disorder. This 
holds however as long the standard definition of spin cur- 
rent is usedi^i 

A basic issue in spin transport theory is the use of 
an appropriate definition of spin current for systems in 
which the spin is not a conserved quantity. The standard 
definition of the spin current operator for spin- 1/2 par- 
ticles is the expectation value of the anticommutator of 
the velocity and spin operators, J s = (j{<t z ,v}), where 
a z is the Pauli spin z— component and v the electron 
velocity operator, respectively. This definition has the 
appealing form that resembles the usual charge-current 
operator, and in addition, for spin-polarized systems, it 
yields to the expected difference between the spin-up and 
spin-down charge currents. However, the conventional 
definition of spin current operator has a caveat, it is not 
conserved in systems with SOI, rendering it incomplete 
in describing a true "spin-current" j 25 i 26 

Recently Shi et aL 27 ^ 2a introduced an alternative def- 
inition of spin current which circumvents the latter is- 
sue. The proposed definition for the effective spin cur- 
rent operator is described by the time derivative of the 
spin displacement operator, which for spin-1/2 carriers 

is simply given by J sz = § d Z** • The alternative def- 
inition adds to the conventional part J s , a spin source 
term (torque dipole density P T ) associated to the elec- 
tron spin precessional motion. In this way the total effec- 
tive spin current density can be defined as J s = J s +P r , 
with J s =Re<I't 1 7 sz iI/. Most important, such effective 
spin current definition straightforwardly satisfies the con- 
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tinuity equation + V • J s = 0, with S z = ^a z V 
describing the spin density. It is also possible to estab- 
lish an Onsager relation between the spin transport co- 
efficients and the mechanical or thermodynamical force- 
driven transport coefficients. In addition, it vanishes for 
localized orbitals, predicting the expected zero spin Hall 
conductivity for insulators^ 

Using the conserved spin-current operator of Ref . [I?! , 
T.-W. Chen et aZ.— predicted that the dc value of the 
spin Hall conductivity for linear Rashba and Dresselhaus 
SOI models has an opposite sign to that obtained using 
the conventional definition. N. Sugimoto et a/.— have 
also studied the SHE and the conditions for nonzero spin 
Hall currents using the conserved current definition for 
both, linear and cubic Rashba models. Subsequently, 
we explored (Ref.[3l|) the spin Hall conductivity for 2D 
electron gases with competing fc-linear Rashba and Dres- 
selhaus SOI in the ac frequency regime. There we have 
highlighted the contrasting results when the conserved 
spin current operator is applied to 2DEGs. Most re- 
cently, T.-W. Chen et a/4 2 , have investigated the spin 
torque and spin Hall currents in generic 2D spin-orbit 
Hamiltonian H so =A(k)cr x — B(k)<j y also using the con- 
served spin-current operator and found that regardless of 
the detailed form of the energy dispersion [i.e. A(k) and 
B(k) coefficients], the conserved static (cu = 0) spin Hall 
conductivity changes its sign with respect to the conven- 
tional spin Hall conductivity. 

In this paper, using the definition for the spin-current 
operator reported by Shi et al^ we explore the behavior 
of the frequency dependent spin (Hall) a s *{u), and of 
the charge conductivity tensor ct^(w) for a 2DHG with 
fc-cubic Rashba and Dresselhaus SOI. We show that the 
optical spectrum of the spin conductivity exhibit remark- 
able changes when this new definition of spin current is 
applied. A rather large response of the spin Hall con- 
ductivity is predicted to arise when using the conserved 
definition of spin current operator owing to a dominant 
contribution of the spin-torque term. In particular, we 
predict that the magnitude and sign of the dynamic spin 
current strongly depends on the electric field frequency 
as well as with the interplay of the Rashba and Dres- 
selhaus spin-orbit coupling strengths. Such behavior is 
similar to that reported in 2DEGsf2i however, unlike the 
later case, which gives vanishing SHE {<jsh{w) = 0) when 
the Rashba and Dresselhaus parameters have the same 
strength, a rather finite spin conductivity (cfsh(^) =h 0) 
is obtained in general for the case of heavy-holes. Use- 
ful analytical expressions for the charge and spin (Hall) 
conductivities as a function of the frequency can be de- 
rived in the common limit e so /sF *C 1, being e so the 
characteristic spin-orbit energy and e f the Fermi energy 
of the heavy-holes. Interestingly, a straightforward con- 
nection between spin and charge conductivities can be 
established in such limit. 

The remaining of the present work is organized as fol- 
lows. In Sec. II the Hamiltonian model for a 2DHG in the 
presence of both, fc-cubic Rashba and Dresselhaus SOI is 



described. The main features of the frequency-dependent 
Kubo formula in linear response and its application to the 
spin-Hall conductivity are presented in Sec. III. We of- 
fer a discussion on the connection between the spin and 
charge conductivities in Sec. IV. Section V is devoted for 
the discussion of the numerical results. We conclude in 
Sec. VI, and finally, in Appendixes A and B, we outline 
the derivation of the spin-current-charge-current corre- 
lation function and the longitudinal spin conductivity, 
respectively. 

II. HAMILTONIAN MODEL 

We are interested to model spin (Hall) transport in 
2DHGs formed in III-V semiconductor quantum wells 
with a relatively strong confinement potential ((fc 2 ) > 
fc 2 ). In such systems, the resulting large splitting of 
heavy hole (HH) and light (LH) hole bands are ex- 
pected to yield fully filled LH states, as reported in re- 
cent experiments^ Hence the relevant contribution to 
the charge and spin conductivity at the Fermi energy 
arises from the HH bands alone. In such case we can de- 
scribe the two-dimensional single particle (spin-3/2) hole 
system by the effective Hamiltonian 

H=^ + H R + H D , (1) 
Am* 

where m* is the effective mass for the HHs, and we 
have included the two most likely dominant SOI terms; 
namely, the fc-cubic Rashba and Dresselhaus spin-orbit 
coupling for HHsj 33 ' 34 In {1} the term Hr denotes the 
Rashba spin orbit coupling which arises from the struc- 
tural inversion asymmetry (SIA) of the hole confining 
potential. It is given by^ 

Hn^-^ia+pl-a-pl), (2) 

with, p± = p x ±ip y , being p x ^ y the components of the 2D 
momentum operator, a± — (a x ± ia y )/2, where a x ^ y are 
the usual spin Pauli matrices, and a specifies the electri- 
cally tunable Rashba SOI parameter for HHs. Values of 
the order of 10 -22 eVcm 3 for a have been found in GaAs- 
AlGaAs quantum well samples with heavy-hole densities 
between 1.8 x 10 10 cm" 2 to 4.2 x 10 10 cm" 2 ^ 

The third term in (fl]), Hp, is the Dresselhaus spin- 
orbit interaction for HHs which can be described bj* 3 ^ 



H D = --^{v+p-p+p- + cr-p + p^p + ). (3) 

This term originates from the bulk-induced inversion 
asymmetry (BIA) of zincblende semiconductor struc- 
tures. The spin-orbit parameter (3 which is fix for a 
given system, can be estimated through the expression, 
P = ilal{p 2 z ){E g + E so )/E so E hh where 7 is the BIA 
(Dresselhaus) parameter, 70 is the Luttinger constant 
(within the spherical approximation, 72 — 73 = 7o)- The 
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FIG. 1: (Color online) Fermi contour kp^(0) showing the an- 
gular anisotropy of the spin-split heavy-hole bands of a 2DHG 
with fc-cubic Rashba and fc-cubic Dresselhaus spin-orbit in- 
teraction. The dotted C r (u>) curve results from the condition 
e+(k) — £-(k) = hu, see text for details. 



energy parameters are as follows, E g is the fundamen- 
tal gap, Ehi and E so are the HH-LH energy gap and 
the split-off-light-hole energy gap, respectively, all at the 
quantum well region. For a GaAs-based quantum well 
of 100 nm of width, (3 is numerically estimated to be of 
the order of 10 -22 eVcm 3 . The Dresselhaus term ([3]) was 
earlier introduced by Bulaev and Loss^l in a study of spin 
relaxation and decoherence in quantum dots in perpen- 
dicular magnetic fields. The importance of the fc-cubic 
SOI related terms in the spin-splitting of the subband 
spectrum of HHs in quantum wells was first noticed by 
Rashba and Sherman^ 

The eigenstates for the full Hamiltonian Eq.JT]) are 
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fie 



1 

i(20-0) 



(4) 



with k = (k x , k v ) = 

ak x —(3k, 



(k cos 9, k sin 9) in polar coordinates, 

tan -1 ^rr-^nr- , A is the area of the system, and 

/i G { + 1- —1} denotes the pseudo-spin ±3/2 HH branch. 
The energy spectrum is cubic in k and it is given by 



h 2 k 2 



2m* 

A(0) = V" 2 +/3 2 - 2a/?sin(2<9). 



(5) 
(6) 



wave numbers kp^(9) which can be written compactly as 



V4u - 3m 2 



LIU 



vM2 



(7) 



being u{9) = 1 - y/l - 2e 2 X{9), with \{6) = A 2 {9)/{a 2 + 

/3 2 ), e = e so /eF and e so — \/ a 2 + (3 2 k F is the charac- 
teristic spin-orbit energy4°-The Fermi energy for vanish- 

h 2 k 2 

ing spin-orbit coupling is estimated as ef = 2m f , with 
Fermi wave number fcf = \J2ixnh and rih is the heavy- 
hole density concentration. 

Now, in the presence of an oscillating electric field of 
frequency u> the allowed direct transitions between the 
HH spin-branches e±(k,9) should satisfy the equation 
e+{k, 9)—E- (k, 9) — ho;, which for a fixed photon energy, 
defines a elliptical curve C r (u>) in the rotated k— space 
described by the equation 

k x {k x + k 2 ) 2 k 2 (k 2 + k 2 ) 2 



1, 



(8) 



with semi-axis of lengths k^(u>) = Hui/2\a — /3\ and 
fc 3 (w) = hui/2\a + j3\ oriented along the principal axes 
(1,1) and (-1,1) of the reference fc-space (Fig. [TJ. The 
symmetry points k a ,b(^>) are given by k a (uj) — fe^_(7r/4) 
and kb(u>) — fci? >+ (37r/4). From this equations two char- 
acteristic energies are determined, 

m*kF 



huj a = 2\a- I3\k 



Huih = 2 1 a ■ 



■p\k% 



1 - 



h 2 

m*kF 



01 

•01 



(9) 
(10) 



additionally, another two absorption frequencies can be 
identified, namely, 



hu± = 2|a=F/3|fc| 



it 



h 2 



(ii) 



which corresponds to transitions between states at the 
points k a (u>) = fci? + (-7r/4) and = /jf,-(37t/4), re- 

spectively. These four characteristic frequencies above 
will play an important role in the overall response of the 
charge and spin (Hall) conductivities as we shall discuss 
later. 



III. SPIN CONDUCTIVITY IN THE 
FREQUENCY DOMAIN 



Notice that, in spite of the obvious difference in the 
spin-orbit Hamiltonian between electrons (linear in k) 
and heavy-holes (cubic in k) in 2D systems, the angu- 
lar anisotropy in the energy spectrum introduced by the 
expression A(9) acquires formally an identical algebraic 
form to that for electrons. 31 ' 39 Nevertheless, the fc-cubic 
spin-splitting of the HH spin-branches, e+(k) — £-(k) = 
2A(9)k 3 , reveals already rather distinct features than the 
case for electrons as discussed below in more detail. 

Following Schliemann and Loss^S for a given Fermi 
energy Ep and vanishing temperature, the two disper- 
sion branches (k) produces two (spin-dependent) Fermi 



Consider a weak and spatially homogeneous electric 
field E(u)y oscillating with frequency to in the plane 
of the 2DHG. Within the linear response theory, the 
frequency-dependent spin conductivity describing a z- 
polarized-spin current flowing in the v = x, y direction 
can be characterized by the Kubo formula 40 




x {J^(t),v y (0)]\^Ar))dt, (12) 
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where e = — |e| is the hole electric charge. Here we have 
assumed noninteracting carriers with /(e^) the Fermi- 
Dirac distribution function (in the limit of zero temper- 
ature) and uj = u> + ir). The vanishing parameter r\ > 
is just an artifact to regularize the integral and guaran- 
tee causality properties of the Kubo formula. 40 However, 
phenomenologically it can be understood as a measure of 
the hole momentum dissipation effects due to impurity 
scattering events, phonons or any other many-body inter- 
actions effects. In this way, ij^ 1 = r will represent here a 
finite life-time of the (spin-orbit coupled) quasi-particles 
due to scattering with holes 4^ 

We remark here that instead of employing the con- 
ventional spin current operator in the Kubo formula, we 
have employed the alternative (conserved) spin current 
operator introduced by Shi et al£± written in the inter- 
action picture, J7^ z (t). The conserved spin current op- 
erator for HHs in the Schrodinger picture is defined by 
J„ z = -|p d dt*^ , where the prefactor 3/2 comes from 
the projection of the total angular momentum of the HH 
states along the growth direction^ To obtain explicitly 
j„ z (t) we first take the ^-component (v = x,y) of the 
current operator with a spin moment polarized along the 
z-axis at time t = (via Heisenberg equation of motion) , 
yielding the operator J° z (0) = J£(0) + JJ(0), where 



*»>-"{'■•-&}■ 



(13) 



is just the conventional spin-current operator definition, 
and 



-TT(r\\ 371 f V vx 

4 m* 



represents the spin-torque contribution to the spin- 
current operator, being V vv * = {y, IV} with T x — 
m*[ap y (3p 2 ~Py) — f3p x p 2 ]/H 4 and {, } denoting the an- 
ticommutator. To get T y simply replace p x — > p y and 
Py — > Px (see Appendix A). Note that in the absence of 
SOI, the operator V v y — 0, and the total spin-current 
operator reduces trivially to the conventional form (fT"3|) . 
The effective spin current operator is then expressed in 
the interaction picture, J* z (t) = e lHt / h J° z (0)e- lHt / h 
which allows us to calculate the spin-current-charge- 
current correlation function of the Kubo formula once 
the single-particle operator of the carrier velocity v y (0) 
is obtained (see Appendix A). After some straightforward 
algebraic manipulations, it is then possible to express the 
frequency dependent spin conductivity as the sum of two 
terms, namely 



(15) 



where the first term to the right, a^ y (Lo), comes from the 
conventional part of spin-current definition and it is given 



by 



3e 



47T 2 7 



2tt 



,(9)- 



cos 9 



dk- 



with g v (6) = (f3 2 - a 2 - 2A(6») 2 )(^ 



(16) 

3 9 + 8 vy sin 9) 

where <5„„/ is the usual Kronecker delta and u>k = 
2Afc 3 /h. The second term in (fTS"]) arises from the spin- 
torque contribution to the net spin current and reads 



6e 



2tt 

dSg v (9) cos6»A 



dk- 



k 10 



(17) 

The ^-integrals in expressions (fT6]) and (fl7|) can be cal- 
culated exactly, however they lead to rather cumbersome 
expressions and shall not be given here. In general, for 
the case a.,/3^0, the 9 integrals can not be performed 
straightforwardly and a numerical integration has to be 
implemented. For a pure Rashba ((3 = 0) or Dresselhaus 
(a = 0) system, the angular dependence of the integrand 
above reduces significantly, leading e.g. to a vanishing 
longitudinal spin conductivity. Although the transverse 
spin (Hall) conductivity yields close analytic form for 
(T xy (u)) and a xy (u), the resulting expressions still some- 
what complicated making almost impossible to retrieve 
any valuable physical insight from them. Thus, before 
going into the numerics, and in order to make progress 
understanding qualitatively the physics here, it is useful 
to consider first the behavior of the spin conductivity in 
the limit e so /ep <C 1. This is a reasonable limit which 
holds typically in 2DHG's in III-V based semiconductor 
heterostructures and for which analytical formulas of the 
frequency-dependent spin conductivity can be derived. 
From Eqs. (fl6|) and (fTTl) . we obtain to leading order in 
EsoIef = e, 



sz,0 



r c,0, 
>'<) 



M+a^H+0( e d ), (18) 



where the dominant contribution from the conventional 
part of the spin Hall [y = x) reads (see Appendix B for 

v = y), 



9e/87r 



n M *M 



YlMl-n^Y' 2 



(19) 



being £ M — 2(a + p(3)kp, while leading contribution of 
the spin-torque can be written in the form 



V) 



tC,0/ 



x-,H 2 
^7^-3^ H ' 



9e/8?r 

where we have defined the auxiliary function, 



(20) 



5 



2,~,2 



(2w+n„y[s : 



2,~,2 



4e 2 



(21) 



It is illustrative to study the behavior of the spin-Hall 
conductivity in the static limit {ui = 0). In the pres- 
ence of weak disorder, it must be analyzed comparing 
the characteristic spin-orbit energy e so with the energy 
scale of the impurity scattering tvq (related to the time- 
relaxation rate r as r\ = r _1 ) by taking into account 
both limiting cases, e so <C tvq and e so 3> tvq. Assuming 
first that the impurity scattering dominates over the spin- 
orbit coupling, we expand eqs. (|19j) and (|20j) in powers 
of e so /hrj <C 1. To lowest order we get 



_c,0 



(0) 



9e/ 



2e 



- D 



(22) 



while crl;°(0) 



, ,, , . . , — 0, being er = akp and Ed = (3k p. 



On the other hand, if the impurity scattering is weak 
compared to the spin-orbit coupling [e so /hr) 3> 1) we 
obtain to second order 



2 1 

"3 + -3 Sgn[a 



9e/8n 

with a nonzero torque part 



2 P 2 ) 



(23) 



9e/87r 



-sgn(a 2 - /3 2 ) 



fcV K 



(24) 



where we have defined the dimensionless parameter 1Z = 
(a 2 + f3 2 ) (3a 2 - (3 2 ) / (a 2 - (3 2 ) 2 . The above expressions for 
c X y(0) reduces, in each limit case, to the known formulas 
for (3 = (only cubic-Rashba) reported in Ref. [4(| ■ 

Moreover, note that, at zero frequency and for ultra- 
clean (77 — > + ) samples, i.e. for Q — ► + i0 + the 
function G(lu) — > 0, which in turn simplify (|20|) to 
= — 2cr^y (0). Thus, the torque-dipole contribu- 
tion makes the total spin-conductivity, a X y ,0 (Q) to change 
its sign with respect to the conventional result. It can be 
shown that this relationship is valid as well for the lon- 
gitudinal spin conductivity, <Jy^(u>) in the same limit. 
From Eqns. ([11]), $Z5§ and ((23$ the static value of the 
spin Hall conductivity for the 2DHG (free of disorder) 
takes the universal form 



-9e/87r for a 2 > (3 2 , 
(r«(0) = { -6e/87r for a = (3, 
-3e/87r for a 2 < (3 2 . 



(25) 



In contrast with the 2DEG case^-^ Eq. (gSJ) pre- 
dicts a non- vanishing dc spin-Hall conductivity (— 6e/87r) 
when the spin-orbit parameters, a and (3 have the same 



strength. There is also an asymmetry of the dc spin Hall 
conductivity depending whether a 2 > /3 2 or a 2 < (3 2 , 
unlike its counterpart for electrons in the clean limit. 
We obtain cr^(0) = -9e/87r for j3 = whiles cr^(0) = 
— 3e/87r for a — 0, in striking difference to the analogous 
for electrons where ±e/87r is obtained, respectively. 



It is known that for a 2DEG in presence of Rashba and 
(linear) Dresselhaus SOI, with spin-orbit parameters a e 
and p e , respectively, the operator S — (a x ± a y )/y2 pro- 
vides an additional conserved quantity when a e = ±/3 e 
since the Rashba and Dresselhaus SOI Hamiltonians for 
electrons commutes with S in such case.— This symmetry 
leads to fc-independent spin-states which in turn yields a 
suppression of the spin Hall effect at a e = ±/3 e . For a 
2DHG with cubic Rashba and Dresselhaus SOI the sce- 
nario is completely different since such symmetry is for- 
bidden for a = ±/3. Indeed, the operator S does not com- 
mute with the 2DHG-Hamiltonian with cubic Rashba 
and Dresselhaus SOI (incidentally this only occurs for 
p y = which yields no spin Hall effect) leading in gen- 
eral to fc-dependent spin states (see Eq. Q). Therefore a 
non-zero response is expected for the heavy-hole spin Hall 
conductivity even for a = ±/3. Alternatively, for elec- 
trons the unitary transformation a x — * — & y , a y — * —o x 
and o z — > — <Jz makes the Rashba and Dresselhaus cou- 
plings to be interchanged, along with a change of sign of 
the spin current, J I — > — J 1 This symmetries explains 
why the spin Hall conductivity only differ by a sign de- 
pending upon the competitions of both strengths, as well 
as its suppression at the symmetry point a e = ±/3 e . How- 
ever, for cubic-Rashba J2]) and cubic-Dressclhaus ([3]) SOI 
there is no a mapping of one Hamiltonian to the other 
(the symmetry above for electrons is broken here), and 
thus, the dc spin Hall conductivity is asymmetric depend- 
ing upon the relative ratio between a and 0. As we will 
see below, the ac spin Hall conductivity for heavy holes 
is also finite for a = ±/3. 



Furthermore, in the ac domain (u> 7^ 0) it can be 
demonstrated that the torque spin Hall conductivity PD|) 
is related to the conventional spin Hall conductivity (fl9|) 



through the relation aZ v (&) = — (2 + ibd / dQ)a% y {Cj) , in 
agreement with Ref.[32j. This relation does not hold 
however for the longitudinal (along the driving field) spin 
conductivity. Moreover a spin current can be generated 
along the driving field direction as long the system is 
submitted to a non-zero Dresselhaus SOI, as has been 
previously reported for electrons^ 
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IV. CONNECTION BETWEEN CHARGE AND 
SPIN CONDUCTIVITIES 

In a 2DHG subject to an oscillating electric field, the 
driven electric current is described by the charge cur- 
rent conductivity tensor T, uy (uj) — 5 uv ctd(u)) +a°y(Lj), in 

which an(uj) = m Z ~ — -r- 2 ^ — , is the dynamic Drude 

2 

conductivity, being a — ^f^r the static Drude con- 
ductivity and af,y(uj) is the contribution due to inter 
spin- split induced transitions. Within the linear-response 
Kubo formalism, following a similar procedure to previ- 
ous section, the spin-orbit induced contribution to the 
charge current conductivity takes the form 



e 2 /2Tth 



2i 



d6h v {0) 



cos t) 



dk- 



k F. - 



(26) 

where h v {9) = {5 VX sm6-5 vy cos<9)(2A 2 +a 2 -/3 2 ) 2 . Note 
that, as it occurs with the dynamic Drude conductivity, 
the model leads to a charge current conductivity which is 
singular for the simultaneously clean (rj — > 0) and zero- 
frequency limit. Although for a,j3 ^ the fc-integrals 
in the above expression are elementary, the subsequent 
integrals in 6 do not have an analytical solution. It is 
then illustrative to examine the limit e < 1 to obtain 
analytical expressions of a^z{ui). It can be shown that 
for a and (3 different from zero, the charge conductivity 
can be related tc 
Appendix B) as 



can be related to eqn. (flU)) for a%^{to) and cr^(w) (see 



T c,0 



9e/87r 
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9e/87r 



(3a 2 + P 2 )k F 



e 2 /2irh 



ch,0 



hhJE} 

2a(3k 6 F 



e 2 /2nh 



frbJEl 



,(27) 



(28) 



being A{u) = (>±/?>)frj0j£ F (2h 2 Cj 2 + Yl^^)' 1 - From the 
equations above it follows that the longitudinal spin con- 
ductivity can be determined through the (transverse) 
charge Hall conductivity induced by the SOI. A similar 
situation occurs for the spin Hall conductivity, which in 
turn can be written in terms of the longitudinal charge 
conductivity. Such connection can be simply interpreted 
as a manifestation of the inverse SHE. It can also be 
viewed as a consequence of the existence of a charge 
imbalance in the sample due to a mean SOI induced 
transverse force in the (anisotropic) fc-cubic Rashba- 
Dresselhaus hole system, as recently predicted to oc- 
cur in such systems by T.-W. Chen et al£i> The latter 
implies the presence of a finite Hall voltage (transverse 
charge conductivity) in the absence of magnetic field and 
as long an electric field is present in the SOI system. 
As the electric field E(u>) is assumed here along the y 
direction, the total spin current is given by J sz (io) — 
[^(w)x+a s y *{uj)y]E(uj). Thus from Eqs. and 
(I27() . the spin conductivity tensor and the spin current 
can be obtained, in principle, through the measure of the 



charge conductivity within the frequency domain alone. 
We believe that this may provide an electrical method to 
detect heavy holes spin currents in presence of fc-cubic 
Rashba and Dresselhaus SOI. 



V. NUMERICAL RESULTS AND DISCUSSION 

We start our discussion of the numerical results by 
considering the isotropic case first, i.e. when only one 
type of SOI, Rashba (f3 — 0) or Dresselhaus (a = 0), 
is present. We considered a 2DHG formed in a GaAs- 
AlGaAs quantum well with a heavy hole effective mass 
of to* = 0.51 m and a moderated sheet hole density of 
rih = 3x 10 11 cm -2 ; here m is the free electron mass. For 
this system, the Rashba parameter has been calculated 46 
to be a — 7.48 x 10~ 23 eVcm 3 , which gives rise to a HH 
spin-splitting at the Fermi energy of A# = 2akp ~ 0.38 
meV. The Fermi wave number for vanishing SOI is es- 
timated from &f = \f2~wn~h . The parameter describing 
the momentum relaxation rate has been chosen such that 
fir) = 0.035 meV, value that would corresponds to sam- 
ples with mobilities [i ~ er/m* ~ 20 m 2 /Vs and relax- 
ation times of r ~ 118 ps. 

In Fig. [2] we plot the numerical integration of the real 
part of the spin-Hall conductivity against the frequency 
of the applied electric field, namely, expressions (fTS")) , (TTBl) 
and (T7]) for the case with only Rashba SOI (/3 = 0). The 
short-dashed (green) line depicts the conventional term 
of spin Hall conductivity, a^ y (uj). It shows a resonance 
like behavior at the energy ftw 
like behavior at hui- ~ 2akp 
would correspond to the minimum and maximum photon 
energy required to induce optical transitions between the 
initial fi = +1 and final /i = —1 pseudo-spin-3/2 split 
branches (see Fig. [T]) . At low frequencies it approaches 
the universal value 9e/87r, while it vanishes for high fre- 
quencies (u> >> Ar/K). Because of the finite value of 
the damping parameter rj, the spectrum has a somewhat 
smoothed shape. 

The torque dipole contribution cr^ y (a;), large-dashed 
(red) lines, shows a rather different behavior. It devel- 
ops a strong resonance at hw ~ Ar as well as near the 
frequencies lu±. At low energies, the torque contribution 
approaches to the value — 2(9e/87r), and as a result, the 
total spin Hall conductivity changes its sign relative to 
the conventional result reaching the value — 9e/87r in the 
static limit. In general, we notice that the frequency- 
dependent torque-dipole response dominates the shape 
of the spectrum of the total spin Hall conductivity and 
consequently, a dramatic change of the overall shape of 
the spectrum <J X y(uj) is observed. Such large effect of the 
inclusion of the torque term in the new definition of the 
spin current has been also reported to occur in 2DEG 
with Rashba and Dresselhaus SOI. 31 The inset of Fig. 2 
shows a comparison of the spin conductivities between 
the case with j3 = ( a = 7.48 x 10 -22 eVcm 3 ), and the 
opposite case, i.e. a = (with (3 = 7.48 x 10 - 22 eVcm 3 ) 



I- ~ 2ak F , and a Fano- 
F _, which from Eq. (fTTj) they 
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FIG. 2: (Color online) Real part of the spin Hall conductivity 
as a function of the photon energy hu for a GaAs-AlGAs 
2DHG system with Rashba SOI only (J3 = 0). The dotted 
(green) curve is obtained using the conventional definition 
of spin-current, whereas the dashed (red) and solid (blue) 
curves are the torque contribution and total spin conductivity, 
respectively, employing the conserved spin-current operator. 
In the inset, the case with j3 = (a = 7.48 x 10~ 23 eVcm 3 ) 
is contrasted with the a = (/3 = 7.48 x 10' 23 eVcm 3 ) case. 
See text for values of remaining parameter used. 



case. Qualitatively, they share essentially the same spec- 
tral features, but with a clear difference in the intensity 
for the central peak. In the dc limit, the effective spin 
Hall conductivity for the system with a ^ and (3 = 
converges to — 9e/8ir, whereas for the case with a = 
and (3 ^ it reaches the constant value — 3e/87r. Such 
asymmetric response is expected due to the formally non- 
equivalent fc-cubic (Rashba and Dresselhaus) Hamiltoni- 
ans, in contrast with the case of electrons, where the 
Rashba and Dresselhaus Hamiltonians can formally be 
mapped into each other, as discussed in section III. 

New interesting features appear when the interplay of 
the Rashba and Dresselhaus SOI is considered. In Fig. [3] 
(a) and (b) we plot the frequency dependent spin Hall and 
longitudinal conductivity for (3 = 0.5a (the remaining pa- 
rameters are as in FigJ^J ■ The results obtained through 
the approximated formula (|18p and the exact numerical 
integration of (|15[) are presented for comparison, showing 
a good agreement, mostly at low frequencies. Here it is 
also evident the remarkable difference between the opti- 
cal spectrum resulting from the use of the standard and 
the conserved spin-current operator definition, J„ ,z . The 
physical origin of the main spectral features can be un- 
derstood in terms of the anisotropic spin splitting caused 
by the simultaneous presence of the Rashba and Dressel- 
haus couplings, in an analogous way to that of the 2DEG 



case 



31,39 




In the limit of vanishing temperature, the sum over 
states in eq. (fT2l is restricted to the region between 
the Fermi contours kp^ + {9) < k < k F _(9), for which 
e+(k,9) < £f < £-(k,8). The number of direct tran- 
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FIG. 3: (Color online) Frequency dependent spin (a) Hall 
and (b) longitudinal conductivity for a 2DHG system with 
finite Rashba and Dresselhaus SOI (/3 = 0.5a). The dotted 
(green) curve shows the result using the standard definition of 
spin current, the case for the conserved spin-current operator 
Jv* in the limit e < 1 is shown in dashed (red) line, and 
the outcome of the exact numerical integration of Eq. (|15[) is 
shown with solid (blue) lines. Here hr\ = 0.035 meV and other 
parameters are as in Fig. [5] We identify four characteristic 
frequencies, two defining the optical absorption edges, lo- and 
oj+, while the other two correspond to the peaks of the spin 
conductivity occuring at ui a and uib- The latter two arises 
due to the symmetry of the spin-split conduction bands in 
k— space at the Fermi level. 



sitions that can take place at the energy Hu) only in- 
volve states with wave vectors that satisfy the equa- 
tion e+(k,6) — E-(k,8) = tkj. As discussed in the in- 
troduction, four distinctive frequencies can be identified 
given the anisotropic k— space available for the optical 
response; namely, uj a , Ub and w±. 

The peaks observed in the optical conductivity cr**(cj) 
of Fig. 3 correspond to energy transitions hjj a ,b involving 
states in the vicinity of the symmetry points k a (uj) = 
/cf.-(tt/4) and kb(to) = fci? j+ (37r/4), respectively (see Eq. 
©)• Similarly, we can see that there are absorption edges 
at energies corresponding to transitions between states at 
the points k a (u>) = kp^ + (Tr/4) and kf,(u>) = fci? j _(37r/4), 
which would correspond to the transition energies Hu)±, 
respectively. For clarity and to guide the eye, such char- 
acteristic energies are indicated in Fig. 3 with dashed ver- 
tical lines. It can also be observed a large separation be- 
tween resonances hwb and hw- unlike the 2DEG case. 31 
Surprisingly, such observation indicates that there is a 
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larger splitting between the Fermi contours kp j+ and fcjr _ 
along the (-1,1) direction for 2DHGs than for 2DEGs. 
The finite value of hrj chosen here, besides of introduc- 
ing an overall smoothing of the spectrum, it also yields 
to a slight shifting of the appearance of the peaks in re- 
lation to these (four) characteristic energies. It is clear 
that the definition of spin-current operator by Shi et al. 27 
yields a drastic different frequency response from that 
predicted by the conventional definition. In addition, as 
it occurs with the pure Rashba (or Dresselhaus) SOI case, 
the torque dipole contribution turns out to be the domi- 
nant term in the spin-Hall conductivity. 

We have also explored the effect induced of varying 
the ratio a/ (3 on the conserved spin conductivity as a 
function of the exciting frequency. In particular, in Fig. 4 
(a) the spin Hall response is shown for the specific values 
of a/ (3 = 0,0.5,1 and 1.5 while fixing the rest of the 
sample parameters as in Fig. 3. Notice that the energy 
separation of the resonance peaks becomes larger as the 
aspect ratio a/j3 is increased. A non-zero ac spin Hall 
conductivity is observed in general for the a — (3 case, 
in contrast to the 2DEG which is zero at all frequencies. 
It is also appreciated that the resonance frequencies w+ 
and u>b disappear for this case. This can be explained 
in terms of the collapse of the equations © and (fTTj) at 
a = (3. The latter is due the overlapping of the spin 
±3/2 dispersion laws, £+(k) = £-(k), at kp± — 7r/4 and 
at Uf± — 37r/4, which occurs precisely at the symmetry 
point a = (3. Such effect is emphasized in Fig. 4(b) where 
a color map of the spin Hall conductivity is plotted as a 
function of a continuous variation of the ratio a/ (3 and 
the exciting frequency. 

In figure 5 we have plotted the (a) longitudinal (v = y) 
and (b) transversal (y = x) spin-orbit contribution to the 
charge conductivity for a 2DHG system with (3 = 0.5a. 
The dashed curve corresponds to the plot with expres- 
sions in the limit e so /ef <C 1, whiles the solid line rep- 
resent those obtained by exact numerical integration of 
Eqn. (f26|) . The results shows a quite good agreement be- 
tween both curves for frequencies smaller than cob- The 
four main frequencies ui + , u> a , ^>b and w_ are shown 
here as vertical lines. Interestingly, we have noticed that 
these distinctive resonance frequencies, can in principle, 
be used to estimate the strengths of the spin-orbit pa- 
rameters a and (3. This is done as follows. From Eqs. © 
and (fTTj) assuming e so /ef <C 1 and a, (3 > we arrive to 
the useful expressions 



aQ(a - (3) + [3Q([3 - a) = 



/3@(a - (3) + a@((3 - a) 



hk. 



8^4 (29) 

Ul + + Ul a — LOb — UJ- 



6m* u + + u a + ujb + UJ- 
(30) 



where &(x) is the usual Heavisidc step function. Assum- 
ing that such characteristic frequencies can be experi- 
mentally identified, either via the spin-Hall conductivity 
or from the charge conductivity tensor measurements, 
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FIG. 4: (Color online) (a) Total spin Hall conductivity for 
different aspect ratios of a/ f3 calculated employing the con- 
served spin-current operator J sz . The parameters are as in 
Fig. 2. Notice that the resonance peaks tend to separate in 
energy as the ratio a//3 is increased, (b) Color contour map 
of the ac spin Hall conductivity showing its behavior with a 
continuous variation of the frequency ui and to the relative 
ratio a//3. 



thus the values for both, the Rashba and Dresselhaus cou- 
pling parameters can simultaneously be estimated from 
the formulas above. First it is needed to retrieve if the 
system has a < /3 or a > f3. This can be done by vary- 
ing a gate voltage to change the heavy-hole concentration 
nh, which consequently changes continously the tunable 
a parameter^ Since (3 does not depend on the carrier 



density, then the plot of the right hand side of Eq. 
should give basically a constant output, which inciden- 
tally will coincide with value of (3 if the 2DHG has a < (3, 
see Fig. 6. Then the value of a is determined through the 
expression in Eq. (|3H)) . On the other hand, if a linear be- 
havior of (|29[) with rih is observed, this is an indicative 
that a > [3 and its value gives an estimate for a for a 
given rih concentration. Once a is determined in this 
regime, the value of the constant parameter (3 can be 
readily obtained from Eq. (|30[) . 
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FIG. 5: (Color online) Frequency dependent (a) longitudi- 
nal and (b) transverse spin orbit contribution to the charge 
conductivity for a 2DHG system with fe-cubic Rashba and 
Dresselhaus SOI {fi = 0.5a). The dashed curves shows the 
result in the limit e < 1 whiles the solid curves corresponds to 
the exact numerical integration of Eq. (|26p . There is a good 
agreement between both cases for frequencies lower than u)b- 
The sample parameters are the same as in Fig [2] The four 
main frequencies uj+, ui a , (±>b and ui- are shown as vertical 
lines to guide eye. 
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FIG. 6: (Color online) Plot of the sum of the frequencies u>+, 
bj a , <Jb, and u>- in units of 8h~ 1 k% as a function of the heavy 
hole concentration, n^. For a > j3 a linear variation with n/, 
is observed, while for the opposite case {(3 > a) an almost uni- 
form behavior is seen in a wide range of hole concentrations. 
The latter responds to the fact that typically a is dependent 
on the heavy holes density, while j3 is not. The dramatic dif- 
ference of the behavior between the cases a > /? and /3 > a as 
the density is varied can in principle be used to experimen- 
tally estimate the value of the coupling parameters. See text 
for details. 



VI. CONCLUSIONS 

In this paper we have examined the spin and charge 
conductivity in the frequency regime for a 2D heavy hole 
gas with /c-cubic Rashba and Dresselhaus SOI employing 
a recently proposed spin current operator. Our results 
shows that the optical spectrum of the spin (Hall) con- 
ductivity changes substantially when the conserved spin 
current operator is used due to a significant contribution 
of the spin torque term. In the dc limit we find that spin 
Hall conductivity for a pure fc-cubic Rashba is 3 times the 
value obtained for the pure fc-cubic Dresselhaus system, 
being a^(0) = —9e/8ir and cr^ z (0) = — 3e/87r, respec- 
tively. Such anisotropy (symmetry breaking) is under- 
stood in terms of the absence of mapping the Rashba <-> 
Dresselhaus Hamiltonians, in contrast with the case for 
electrons. It is argued that such asymmetry it is also 
responsible for the non-vanishing dc spin Hall conduc- 
tivity (<r^y(0) — — 6e/87r) when the spin-orbit Hamilto- 
nians, Rashba and Dresselhaus have the same coupling 
strength. In the ac limit, it is shown that the magni- 
tude and the sign of the dynamic spin current is rather 
sensitive to the frequency and to the spin-orbit (a and 
j3) coupling strengths. The latter suggests the possibility 
of the optical manipulation of spin currents in addition 
to the control obtained through external bias. The dra- 
matic differences of the spin Hall conductivity response 
predicted here with the use of the new conserved spin- 
current operator in relation to the calculated with its 
conventional counterpart, it is clearly a call for the ex- 
perimentalist to perform measurements in the frequency 
domain and validate (or discard) its applicability in de- 
scribing spin-transport. In addition, we have shown that 
the angular anisotropy of the spin-splitting energy in- 
duced by the interplay between the Rashba and Dressel- 
haus couplings gives rise to four characteristic resonance 
frequencies, which in principle, can be employed to esti- 
mate the strengths of the spin-orbit parameters through 
optical spectroscopy and/or transport measurements. 
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APPENDIX A: 
SPIN-CURRENT-CHARGE-CURRENT 
CORRELATION FUNCTION 

In this appendix we briefly outline the derivation of the 
expectation value for spin-current-charge-current correla- 
tion function {^(r)\[J^ z {t),v y {Q)]\^{v)) = ^(k,t) 
appearing in the Kubo formula Eq. (fT2")l . We begin by 
writing the carrier velocity operator via the Heisenberg 
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equation of motion v(0) — i[H,r]. For the v 
poncnt, the velocity operator reads, 



y com- whereas the torque spin-current operator yields 



Vy(0) = — 



Py 



dT x 



dT 



^± + a y ^L , (Al) 



dp. 



dp 



where we have defined, 

T x = ^r- [ap y (3pl -p 2 y )~ PPxP 2 



m 



(A2) 



To get T y simply replace p x — > p y and p y — > in (|A2|) . 
It follows trivially that in the absence of SOI, the oper- 
ator T XtV — 0, and the total velocity operator reduces to 
the dispersionless ioimp y /m* . Next we proceed to calcu- 
late the operators x(t), a x (t) and a y {t) in the interaction 
picture, 6{t) = e lH ^ h 6(Q)e- tHt ^ h . After some algebraic 
manipulations, the position operator reads, 



t f 9T~ OL v \ 



+ 



hp y ((3 2 - a 2 



dpx 
2A 2 ) 



2p 2 



v dp x . 
JC(t), 



(A3) 



whiles the spin Pauli operators in the interaction picture 
take the form, 

<r x {t) = <t x {0) + (i) , a y (i) =a v {0)- Yj v (t) , ( A4) 

here we have introduced the time-dependent operators 



£ x (t) = m*T y C(t) 
E y (t) = m*T x C(t) 



(A5) 
(A6) 



huj, 



[ujpt — sm(ujpt)} lZ xy — 2sin 2 (uj p t/2) <j z 



■| sin(w p i)er 2 - — — 7£ ay sin 2 (w p t/2)}, 



(A7) 



in which 
K{t) = 
L{t) = 

lUJJp v. IIWp 

with lZ xy = m*(a x T y — a y T x ) and lu p — 2Ap 3 /ft 4 . Note 
that a Xl a y , p, p x and p y are all given in the Schrodinger 
picture. The expressions above are needed in the calcu- 
lation of the commutator [J x ' z (t),v y (0)], in which the 
effective spin-current operator is written in the interac- 
tion picture, i.e. J^' z {t) = e iHt ' h J^ z {Q)e' lHt ' h . At 
t = we have J°' z {0) = J°(0) + JJ(0), as given explic- 
itly by Eqs. (fT3)) and (I14[) . Then the commutator of 
the velocity operator with the conventional spin current 
reads, 



iHt/h 



l. to j 



iHt/h 



3fi p K 



o-*(t),- 



5T* 



,«v(0) 
<9T 



a 

(A8) 



3ft 
3ft 



e iHt/fi , 


'XX j 


e «ft/ft < 


1 

v TO* J 



-iHt/h 



,0,(0) 



y I -zHt/h t . 



(A9) 



with = {f, IV}. After some algebraic manipula- 

tions we arrive to, 



3ft 



3 2 



3ft 2^2, 
2 m* 



(A10) 



with £ = i(0) and the remainder operators are 



Kf 1 dT x 
Ai = — -t: — cr* , 



(A11) 



TO* 

Ml = — COs(cJpt) CTa; , A^2 = -—Fx COs(cJpt) f7 y , 

(A12) 

sin(uj p t) a z . (A1>!) 



M 3 = ^ sm(uj p t) a z . 



Adding up the expectation value of (|A8[) and (|A10[) 
with the aid of (U]) , we finally arrive to a explicit formula 
for the spin-current-charge-current function ^(k, i), 
which is given by 



, . (k 2 x -k 2 v )k 2 (2a 2 + A 2 ) , 
J>(k,t) = 6^ y - — cos{u k t) 



-3i/i 



Hk 2 x k {I3 2 -a 2 -2A 



2\ , 



cos(cjfet) - cjfc t sin(wfct) , 
(A14) 



with u>k = 2Ak 3 /H and /Lt = ± the pseudo-spin ±3/2. 
The result of Eq. (|A14|) together with the factor e iQt in 
Eq. (p~2|) leads expressions exactly integrable in time for 
the conventional (HHJ) and torque (fT7j) contribution to the 
effective spin Hall effect. 



APPENDIX B: LONGITUDINAL SPIN 
CURRENT 

In this appendix we write down useful analytical ex- 
pressions for the ac and dc longitudinal spin conductivity, 
a y z (ui). Expanding Eqs. (fT6)) and (fT7|) to leading order 
in E so /eF by assuming e so /ef << 1 for a and f3 differ- 
ent from zero, the conventional part of the longitudinal 
(v = y) spin conductivity takes the form 
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-9e/8n 



h 2 Q 2 - 4e 2 so 



2h 



2,~,2 



8a(3k 6 F 



The spin-torque contribution reads 

9e/87r "' 9e/8^ 3 1 
with the function T~L{ui) is defined as 



(Bl) 



(B2) 



Clearly, for the case of identical SOI coupling strengths 
(a = (3) the function H(uj) vanishes at all frequencies 
and crT^ jbj) = —Aa^{uj). It can be shown that using 
Eqs. (|Blj) . (|B2)) and ffe]) within the dc limit and for ultra- 
clean {jj — > + ) samples, the static value of the longitu- 
dinal spin conductivity will reduce to 



i(0) = { 



for 

07T a 



3e a 
8tt /3 



a 2 > (3\ 
for a = (3, 
for a 2 < /3 2 . 



(B4) 



/3 2 ) 



1 - 



h 2 Q 2 - 4e 2 



4a/3 
4a/» 2 [EU, 



h 2 Co 2 Y/ 2 



2h 2 u 2 ]k% 



which depends on the ratio a/f3, as it occurs for a 
(B3) 2DEG4 2 - 
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